Draft version December 1, 2010 

Preprint typeset using I4TgX style emulateapj v. 1 1/10/09 



MASS AND DENSITY OF B-TYPE ASTEROID (702) ALAUDA 

P. Rojo 

Departamento de Astronomia 
Universidad de Chile 
Santiago, Chile 



o 

(N 

> 
O 

o 

m 

6 



> 
in 



J. L. Margot 

Departments of Earth and Space Sciences, Physics and Astronomy 
University of California, Los Angeles 
Los Angeles, CA 90095 
Draft version December 1, 2010 

ABSTRACT 

Observations with the adaptive optics system on the Very Large Telescope reveal that outer main belt asteroid 
(702) Alauda has a small satellite with primary to secondary diameter ratio of ^56. The secondary revolves 
around the primary in 4.9143 ± 0.007 days at a distance of 1227 ± 24 km, yielding a total system mass of 
(6.057 ± 0.36) x 10 18 kg. Combined with an IRAS size measurement, our data yield a bulk density for this 
B-type asteroid of 1570 ± 500 kg m~ 3 . 

Subject headings: Astrometry; Ephemerides; Instrumentation: adaptive optics; Minor planets, asteroids: indi- 
vidual: 702 Alauda; Minor planets, asteroids: general; Techniques: high angular resolution 



1. INTRODUCTION 

The discovery of solar system binaries has a considerable 
impact on key problems in planetary science, partly because 
the binaries allow direct measurements of fundamental phys- 
ical and chemical properties that are otherwise only obtain- 
able with spacecraft. Characteristics of the mutual orbit pro- 
vide crucial information about asteroid bulk properties and 
internal structure, such as mass, density, porosity, and me- 
chanical strength. These measurements are used to estab- 
lish links between asteroids and meteorites, and to under- 
stand the geological context of meteorites that can be stud- 
ied in great detail in the laboratory. Density measurements 
can also constrain the proportion of ice to rock in distant 
minor planets, which is a strong indic ator of the chemical 
environment at the time of form ation (iLunine and Tittemord 
119931: iJohnson and Luninell2005l) . The current proportion of 
binary systems and their configuration are strong tracers of 
the collisional and dynamical evolution of small bodies in 
the solar system. Binaries therefore illuminate the condi- 
tions in the solar nebula and early solar system, and they help 
us refine our understanding of planet formation. The very 
high science pri ority of characterizing binary systems has 
been r e viewed bvlWeidenschilling et"aT[ ( 1989); M erline et al.l 
(2002); Margotl (120021) : iBurnsI d2002l) : iNoll et all (120081) . 

The third largest asteroid and former planet (2) Pallas 
(a=2.772 AU, e=0.231, i=34.8°) belongs to the taxonomic 
class B, which is characterized by a linear featureless spec- 
trum with bluish t o neutral slope i n the wavelength range 
(0.435-0.925 ^m) dBus and Binzell 120021) . In the case of 
Pallas, a broad absorption band in the 3 j jm region is de- 
tected (lLarson et al.ll983l:lJones et al.ll99 0) suggestive of hy- 
drothermal alt eration and remin iscent of carbonaceous chon- 
drite material (ISato et al.lll997l) . The negative spec t ral slo pe 
of B-type asteroids is attributed by Yan g and Jewittl (120101) to 
a broad absorption band near 1.0 /im. Because B-type aster- 
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oids may represent relatively primi tive material with similari- 
ties to carbonaceous meteorit es (|Larson et alJl983HSato et ail 
Il997t I Yang and Jewittl 120101) . their characterization can fur- 
ther our understanding of the primordial building blocks of 
planets. Unfortunately, our knowledge of the density of B- 
type asteroids is poor, in part because few measurements 
exist, and in part because mass estimates for (2) Pallas are 
widely scattered as they rely on indirect techniqu es (gravita- 
tional perturbations on other asteroids or Mars). Carry et al. 
(2010) reported a density of 3400 ± 900 kg m" 3 based 
on a mass esti mate derived fr om pertubations to the orbits 
of asteroids dMichalakl 120001) with 22% un certainties and 
their v olume estimate with 4% uncertainties. Schmi dTet al.l 
(2009) reported a density of 2400 ± 250 kg m -3 based on 
a mas s estimate derived f rom perturbations to the orbit of 
Mars dKonopliv et al. 2006) with 2.7% uncertainties and their 
volume estimate with 10% uncertainties. However their re- 
port lists mass estimates for Pallas that vary by 20%, and 
a more e xtensive surv ey indicates mass estimates th at vary 
by 60% (iHiltonl 120021). Us ing the mass estim ates of iHilton l 
(1999) and Michalak (120001) and the volume of dSchmidt et al.1 
2009) yields densities for Pallas in the range of 3700 ± 390 
kg m -3 and 2800 ± 670 kg m~ 3 , respecti vely. Here we report 
on the discovery dRoio an d Margot 2007) and tracking of the 
first satellite to a B-type asteroid, (702) Alauda (a=3.191 AU, 
e=0.022, i=20.6°), which allows us to provide a direct mass 
measurement and improved density estimates for this class of 
asteroids. 

(702) Alauda has been identified as the largest member of 
a highly i nclined dynamical family in the outer m ain belt of 
asteroids dFoglia and Masill2004l; lGil-Huttonll2006l) . It is rea- 
sonable to assume that the creation of this binary system dates 
bac k to the impact responsible for the f ormation of the fam- 
ily dMichel et al.ll200lUDurda et al.ll2004l) . Because tidal evo- 
lution depends on the formation age of the binary and the me- 
chanical properties of the components, our observations can 
also be used to provide a relationship between formation age 
and mechanical strength, specifically shear modulus and tidal 
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Figure 1. H-band image of the (702) Alauda binary system. The companion 
is clearly seen in the lower-right quadrant. 

dissipation factor (Tayl or and Mar got 2010). 

2. OBSERVATIONS AND DATA REDUCTION 

Using a highly optimized observing strategy, we were able 
to observe a total of 1 1 1 unique main belt and Trojan aster- 
oids during three nights in visitor mode at the Very Large 
Telescope ( VLT). We used the adaptive optics (AO ) instru- 
ment NaCo dLenzen et alj|2003t iRousset et alJl2003l) that has 
an 1024x1024 pixel InSb array detector, which is fed by an 
adaptive optics system consisting on Shack-Hartmann sensor 
with a 14x14 lenslet array. Information from the wavefront 
sensor is used to control both a tip-tilt mirror and a deformable 
mirror. We used NaCo in its S 1 3 mode, resulting in a pixel 
scale of 13.72 mas. 

The observing sequence for each target consisted of four 
10-second integrations in each of four different offset posi- 
tions. Each frame was then flat-fielded, and subtracted in pairs 
to remove the sky, before alignment and addition to form a 
composite image. We also obtained a bad-pixel mask from 
outlier pixels present in the averaged flat-field, and used the 
mask to ignore those pixels in each of the 10-second frame. 

The companion to (702) Alauda was identified in a 40- 
second composite image in the Ks filter (A c = 2.18/im, 
FWHM = 0.35/im). We eliminated the possibility of one 
type of AO artifact (a non-rotating waffle mode) by immedi- 
ately performing an observing sequence with the field of view 
rotated 50 degrees. The non-sidereal tracking of the telescope 
discarded fixed background sources and most moving back- 
ground sources. Follow-up observations that night and the 
next fully confirmed the binarity of the system. 

The binary system was observed at 6 epochs during the last 
2 nights of the observing run with different filters. The H-filter 
(A c = 1.66/im, FWHM = 0.33/im) was found to give the 
best compromise between SNR and AO correction (Fig. [TJ. 
In order to secure additional astrometric measurements, we 
submitted a Director's Discretionary Time (DDT) proposal 
to follow-up the system in queue mode, using the same ob- 
serving instrument and strategy. Table Q] summarizes all the 
observations of the binary system. 

We attempted to detect the companion in our composite 
images, and also in images where an azimuthally averaged 
profile was subtracted from the primary source. We obtained 
ten secure detections. The secondary may be undetectable in 
other images due to proximity to the primary, atmospheric 



seeing, lunar phase, or a combination of factors. We per- 
formed aperture photometry of the primary on the composite 
images, and of the secondary on the profile-subtracted im- 
ages, removing a slanted-plane sky background where neces- 
sary. We then combined the primary and secondary photom- 
etry to measure primary-to-secondary flux ratios and relative 
positions. 

Our orbit determination software uses the separation and 
position angle of the secondary with respect to the primary 
to solve for the orbital parameters in the two-body problem. 
The mutual orbital plane orientation is assumed constant, but 
heliocentric motions and light-time corrections are taken into 
account. Astrometric positions and their errors are specified 
at the mid-time of the exposure sequence. Starting from thou- 
sands of initial conditions, the software adjusts for seven pa- 
rameters (equivalent to six orbital elements plus the mass of 
the system) in the nonlinear ordinary least squares problem 
with a Levenberg-Marquardt technique. The covariance ma- 
trix and post-fit residuals are computed and inspected. Bi- 
nary orbits have been computed wi t h this algorithm f or near- 
Earth asteroids ( Margot et al. 2002; lOstro et al.ll2006l). main - 
belt asteroids (iMargpt and Brownll2003HMerhne et al.l [2002). 
Kuiper belt objects dNoll et alJl2008HPetit et al.ll2008l) . dwarf 
planets, and binary stars. 

3. RESULTS 

Primary-to-secondary flux ratios from the photometric 
analysis are shown in Fig. [2] Some, but not all, of the ob- 
served variations may be related to the rotation of the primary 
with period 8.35 39 ± 0.0007 hrs and lightcurve amplitude of 
0.09 ± 0.02 mag (lBenishekl2008l) or to the rotation of the sec- 
ondary, which we expect to be spin-locked (tidal despinning 
timescale ~10 6 years). Considering the theoretical Poisson 
noise (< 0.05mag), the differences in the same-night ratio 
measurements for the J and H filters (0.21 and 0.27 magni- 
tudes, respectively) must be attributed to unknown system- 
atics. The average of the primary-to-secondary flux ratio is 
3188 ± 1444 (magnitude difference of 8.8tg;y), where the 
quoted uncertainty is the standard deviation of our measure- 
ments. This ratio is one of the largest for solar system bina- 
ries, although a satellite of (41) Dap hne is reported to be 10 
magnitudes fainter than the primary dConrad et al.H2008l) . If 
we assume that the primary and secondary have similar albe- 
does, the primary-to-secondary diameter ratio is ~56, and the 
mass of the secondary represents an insignificant fraction of 
the total mass budget. 

Relative positions of the components at the epochs of our 
measurements are shown in Table |2] Since we only have one 
measurement per epoch, as an upper limit we assigned center- 
ing uncertainties of 0.013 arcseconds in x and y, equivalent to 
the plate scale of the detector. 

Our 7-parameter orbital fits (Fig. 0) used 16 measurements 
at 8 epochs. The best-fit reduced chi-square value was 0.665, 
indicating that our uncertainties were assigned somewhat con- 
servatively but not overly so. Values of the fitted parameters 
and their formal standard errors are shown in Table [3] We 
measure a system mass of (6.057 ± 0.36) x 10 18 kg (roughly 
10~ 6 Earth masses), corresponding to a fractional precision of 
6%, considerably better than results from indirect techniques. 

The IRAS Minor Planet Survey (IMPS) lists the diameter 
of (702) Alauda as 194.73 ± 3.2 km. Combining our mass 
measurement with this size estimate, we compute a density 
with formal uncertainties of 1570 ± 120 kg m~ 3 . Realistic 
uncertainties need to take into account the fact that there may 
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UT Date and time 


MJD 


Filter 


RA (deg) 


DEC (deg) 


A (AU) 


Detection 


2007 07 26 06:24 


54307.2664 


Ks 


316.7047 


-3.7504 


2.1555 


1 


2007 07 26 07:26 


54307.3098 


H 


316.6957 


-3.7485 


2.1554 


1 


2007 07 26 07:35 


54307.3159 


J 


316.6942 


-3.7482 


2.1553 


1 


2007 07 26 07:46 


54307.3235 


H 


316.6927 


-3.7479 


2.1553 


1 


2007 07 27 03:44 


54308.1553 


H 


316.5163 


-3.7124 


2.1528 


1 


2007 07 27 05:58 


54308.2486 


H 


316.4965 


-3.7085 


2.1526 


1 


2007 08 18 04:07 


54330.1716 


H 


311.7363 


-3.2049 


2.1565 


1 


2007 08 29 04:12 


54341.1750 


H 


309.6719 


-3.2051 


2.2082 





2007 09 09 01:45 


54352.0731 


H 


308.1241 


-3.2997 


2.2884 


1 


2007 09 10 02:45 


54353.1147 


H 


308.0073 


-3.3115 


2.2974 





2007 09 11 01:39 


54354.0685 


H 


307.9053 


-3.3225 


2.3059 


1 


2007 09 12 04:09 


54355.1727 


H 


307.7935 


-3.3355 


2.3159 





2007 09 13 02:12 


54356.0920 


H 


307.7056 


-3.3464 


2.3245 





2007 09 14 01:46 


54357.0738 


H 


307.6167 


-3.3581 


2.3338 


1 


2007 09 15 01:10 


54358.0488 


H 


307.5338 


-3.3698 


2.3432 





2007 09 16 02:18 


54359.0959 


H 


307.4507 


-3.3823 


2.3535 





2007 09 20 01:06 


54363.0457 


H 


307.1929 


-3.4286 


2.3940 





2007 09 21 00:54 


54364.0377 


H 


307.1420 


-3.4398 


2.4046 





2007 09 28 01:56 


54371.0804 


H 


306.9412 


-3.5097 


2.4841 





2007 09 30 03:32 


54373.1471 


H 


306.9351 


-3.5261 


2.5086 






Table 1 

Observational circumstances. Times listed refer to the start of the third exposure, roughly the middle of the four-image exposure sequence. Observation epochs 
are given in Universal Time (UT) and Modified Julian Date (MJD) formats. The Right Ascension (RA), Declination (DEC) and geocentric distance (A) are used 
in orbital fits. Observations above and below the line were obtained in visitor and queue mode, respectively. 
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Figure 2. Primary-to-secondary magnitude differences for positive detec- 
tions. Triangle, squares, and circles indicate the measurements using the J, 
H, and Ks filters, respectively. The error bars indicate the Poisson uncertainty. 



be a bias in the asteroid size and that the shape of the asteroid 
may differ substantially from a sphere. The IMPS diameter 
measurements have been compared wi th occupation result s 
and are believed to be accurate to 10% (Ted esco et al.l 2002). 
Therefore we suggest a density estimate with more realistic 
errors of 1570 ± 500 kg m~ 3 . 

The confounding effects of porosity must be kept in 
mind in interpreting d ensity measurements (Britt et aTll2002t 
iMcKinnon et ai]|2008h . In Fig. H we show the range of grain 
densities that are compatible with our bulk density for sev- 
eral values of porosity. Because the pressure at the center of 
(702) Alauda is moderate (~3.3 MPa or ^33 atmospheres), 
this body could in principle sustain fairly large porosities. At 
near-zero porosities a mix of anhydrous silicates and water ice 
would require a fair amount of ice (roughly 3/4 by volume or 
4/9 by mass) to match the density constraint. At high porosi- 
ties the asteroid could be entirely devoid of ice. Progress 
in fully understanding the manifestation of porosity in small 
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Figure 3. Observed secondary-to-primary separations and position angles, 
with error bars (crosses), and best-fit orbit (solid lines). The exterior ellipse 
shows the projection of the orbit on the plane of the sky (62 degree inclina- 
tion) on MJD 54307; it is representative of the first four astrometry data points 
at position angles 21, 25, 58, and 60 degrees. The interior ellipse shows the 
projection of the best-fit orbit on the plane of the sky (69 degree inclina- 
tion) fifty days later; it is more representative of the last three astrometry data 
points at position angles 46, 49 and 213 degrees. 



bodies will require a large sample of reliable sizes and densi- 
ties, or in-situ seismic experiments. 
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UT Date and time 


MJD 


Sep (asec) 


PA (deg) 


Res. (North) 


Res. (East) 


2007 07 26 06:24 


54307.2663 


0.584 


21.3 


0.975 


0.588 


2007 07 26 07:46 


54307.3235 


0.629 


24.9 


-0.081 


-0.167 


2007 07 27 03:44 


54308.1553 


0.751 


58.0 


0.942 


-0.446 


2007 07 27 05:58 


54308.2486 


0.726 


60.5 


-0.063 


0.014 


2007 08 18 04:07 


54330.1716 


0.748 


232.9 


-0.418 


-1.165 


2007 09 09 01:45 


54352.0731 


0.751 


46.4 


-0.924 


-0.554 


2007 09 11 01:39 


54354.0685 


0.651 


213.4 


0.530 


0.357 


2007 09 14 01:46 


54357.0738 


0.724 


48.7 


-0.640 


-0.230 


Table 2 



Astrometry used in our orbital solution, showing the positions of the secondary relative to the primary, expressed as separations (Sep) in arcseconds and position 
angles (PA) in degrees East of North. Epochs are given as in Table[T] Also shown are residuals from our best-fit solution in the format (0-C)/o\ 
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Parameter 




Estimfttc 


^formal 


Orbital period [days] 


P 


4.9143 


0.007 


Semi-major axis [km] 


a 


1226.5 


24 


Eccentricity 


e 


0.003 


0.02 


Longitude of ascending node [deg] 


fi 


280.9 


3 


Inclination [deg] 


i 


46.5 


3 


Argument of pericenter [deg] 




169.074 


180 


Epoch of pericenter passage [MJD] 


T 


54308.648 


4 



Table 3 

Results of orbital fits to our astrometry, giving parameter estimates and 
formal standard errors. Angles are provided in the equatorial frame of epoch 
J2000. The orbit normal is at right ascension 190.9° and declination 43.5°, 
corresponding to ecliptic longitude 168.3° and latitude 43.3°. The 
orientation of pericenter and time of pericenter passage are poorly 
constrained, as expected for such a low eccentricity orbit. Nevertheless we 
provide guard digits in those quantities to allow computation of the relative 
positions of the components without loss of precision. 
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Figure 4. Range of porosities and grain densities admitted by our density 
measurement (our nominal solution is the bold solid line, with density uncer- 
tainties captured by the solid lines above and below). The densities of water 
ice and anhydrous silicates are shown as dashed lines. 



4. CONCLUSIONS 

We report the discovery of a satellite to (702) Alauda, a 
member of an asteroid collisional family. The satellite is 
small, with a ratio of primary to secondary radii of about 
50. The secondary revolves around the primary in 4.91 days 
at a distance of 1,230 km, corresponding to 12.6 primary 
radii. Due to tidal interactions, its spin period is expected to 
synchronize to its orbital period on million year timescales. 
The secondary likely formed by a sub-catastrophic colli- 
sion o n the primary, and subsequently tidall y evolved out- 
ward (|MEnniitiI3llQQl iDurSiliLlEQOl). Because the 
orbit is quasi-circular it is likely t hat the secondary is me- 
chanically weaker than th e primary ( Margot and Brownl2003l : 
iTavlor and Margotl 12010), as might be expected if the body 
formed by the re-accumulation of impact ejecta. Our mass 
measurement with 6% uncertainties provides an important 
new constraint on the composition and internal structure of 
B-type asteroids. Combined with an IRAS size measurement, 
the mass yields a bulk density of 1570 ± 500 kg m~ 3 . This 
density admits an ice-to-rock mass ratio in the range 0-80%, 
depending on the porosity of the asteroid. Our density deter- 
mination is significantly lower than previous estimates for the 
B-type asteroid (2) Pallas, which may be related to different 
porosities. 
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